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Using the adaptive time-dependent density matrix renormalization group method, we numerically 
study the spin dynamics and transport in one-dimensional spin-1/2 systems at zero temperature. 
Instead of computing transport coefficients from linear response theory, we study the real-time 
evolution of the magnetization starting from spatially inhomogeneous initial states. In particular, 
we are able to analyze systems far away from equilibrium with this set-up. By computing the time- 
dependence of the variance of the magnetization, we can distinguish diffusive from ballistic regimes, 
depending on model parameters. For the example of the anisotropic spin-1/2 chain and at half 
filling, we find the expected ballistic behavior in the easy-plane phase, while in the massive regime 
the dynamics of the magnetization is diffusive. Our approach aUows us to tune the deviation of the 
initial state from the ground state and the qualitative behavior of the dynamics turns out to be valid 
even for highly perturbed initial states in the case of easy-plane exchange anisotropies. We further 
cover two examples of nonintegrable models, the frustrated chain and the two-leg spin ladder, and 
we encounter diffusive transport in all massive phases. In the former system, our results indicate 
ballistic behavior in the critical phase. We propose that the study of the time-dependence of the 
spatial variance of particle densities could be instrumental in the characterization of the expansion 
of ultracold atoms in optical lattices as well. 



I. INTRODUCTION 

Transport in low-dimensional strongly correlated sys- 
tems continues to excite theoretical and experimental 
physicists alike. For theorists, transport problems pose 
a formidable challenge, as established tools to work out 
ground-state properties of strongly correlated systems do 
not always provide an adequate description of transport 
as well (see Refs. [H andU and references therein), which, 
in particular, pertains to systems driven out of equilib- 
rium. 

Within linear-response theory, one often distinguishes 
between ballistic and diffusive transport by invoking the 
notion of the Drude weight)^ i.e., the prefactor of a 
delta-function at zero-frequency in the frequency depen- 
dent transport coefficient. A finite Drude weight defines 
ballistic transport, while in the case of a vanishing Drude 
weight, the zero- frequency limit of the conductivity's reg- 
ular part determines the long-time behavior. 

Significant theoretical attention has been devoted to 
one-dimensional spin systems (see Refs. [l|0l5| for a re- 
view). Open theoretical questions include, for instance, 
the finite temperature transport of the anisotropic spin- 
1/2 chain with nearest-neighbor interactions (the XX Z 
chain), with an unsettled debate on whether finite tem- 
perature transport in the Heisenberg chain is ballis- 
tic or not^iii^iSiiiSiiiii^ii^ as wcU as on the actual tem- 



perature dependence of the Drude weight i^ i^^'^^ For 
nonintegrable models and finite temperatures, one ex- 
pects diffusive transport on general grounds ^^li^ii^ 
and numerical studies have widely confirmed this pic- 
ture in the high-temperature limit and massive phases 
of spin models A similar scenario has 
emerged for thermal transport Yet, the issue 
of (quasi)ballistic transport in gapless phases of, e.g., 
the frustrated chain fiSiii^i^i^ at low temperatures is still 
under scrutiny. Moreover, the possibility of anomalous 
transport due to a diverging coefficient of the dc conduc- 
tivity has been emphasized 32i22ii^ 

Much less is known about the non-linear transport at 
large external driving forces, or more generally, nonequi- 
librium properties. A recent study using a quantum 
master-equation approach has addressed the spin trans- 
port in the antiferromagnetic phase of the XX Z chaini^ 
The time-evolution of magnetization profiles in analyti- 
cally exactly solvable models has been the case of interest 
in Refs. H^illliii. 

Besides the fundamental interest in understanding 
large-bias and out-of-equilibrium phenomena, research 
into transport properties of low-dimensional spin systems 
is strongly motivated by exciting experimental results on 
large thermal conductivities in spin ladder and chain ma- 
terials (see, e.g., Refs. [3ll32i for a review). Most evidently 
in spin ladder materials such as (La,Sr,Ca)i4Cu2404i, 
such large thermal conductivities have been attributed 
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to magnetic excitations i^'^'^ In a more recent experiment 
on La5Ca9Cu24 041)2^ heat dynamics has been probed by 
time-of-flight measurements. In this set-up, the surface 
of a sample is covered with a thin fluorescent layer. After 
shining on that surface with a laser, one can then follow 
the propagation of heat in the surface by thermal imag- 
ing at different times. A pronounced difference is seen 
comparing a surface that contains ladders to one that is 
perpendicular to the ladder direction. In the former case, 
heat diffuses predominantly along the ladder direction, 
while little dynamics is seen in the latter case. These 
results support the notion of anisotropic heat transport 
in this material;^ due to the contribution of magnetic 
excitations. 

Besides thermal transport measurements, questions 
of diffusive versus ballistic transport have been experi- 
mentally probed in nuclear magnetic resonance^^*^ (see 
Ref. m for related theoretical work) as well as muon spin 
resonance experiments.— In a more recent development, 
transport properties of low-dimensional ultra-cold atom 
gases have gained attention as well, with experiments 
focusing on the detection of Anderson localization.—!^ 
Interacting two-component Bose gases in optical lat- 
tices have been suggested to potentially realize spin-1/2 
Hamiltonians.— 1^ 

As far as numerical approaches are concerned, on 
the technical side, full exact diagonalization (ED) stud- 
ies are restricted to system sizes of about 24 sites in 
the case of a spin-1/2 chain, while in Quantum Monte 
Carlo simulations, the calculation of frequency depen- 
dent properties of more than two-point correlation func- 
tions remains difficult (see, e.g., Ref. [H and references 
therein). The density matrix renormalization group 
(DMRG) mcthod^i^i^ has most successfully been ap- 
plied to zero-temperature phenomena. With the ad- 
vent of the adaptive time-dependent DMRG method 
(tDMRG),— 1^1^ the study of non-equilibrium and large 
bias phenomena has become possible. While the meth- 
ods mentioned so far are designed for pure bulk systems, 
quantum-master approaches that account for dissipation 
by incorporating baths have been used as well in the 
study of transport in ID spin systemsi^'^'Si^'^ Such 
methods are not constrained to the linear response either 
and circumvent the use of Kubo formulae, as they mea- 
sure gradients and current expectation values directly. 

In this work we will introduce and exploit an alter- 
native approach based on zero-temperature tDMRG cal- 
culations. Instead of analyzing currents and their cor- 
relation functions directly, we study the magnetization 
dynamics after preparing the system in inhoniogeneous 
initial states. For instance, we subject the system to an 
external magnetic field of Gaussian shape and, after re- 
leasing the confining field, follow the time evolution of 
the magnetization. Computing the variance of the mag- 
netization allows us to distinguish ballistic from diffu- 
sive regimes, depending on model parameters: we con- 
sider the dynamics to be ballistic if the variance grows 
quadratically in time, which is the behavior of noninter- 



acting particles, while a diffusive behavior manifests itself 
in a linear increase of the magnetization's variance. Our 
approach has the advantage that we can characterize dif- 
fusion by following the time-evolution of a local quantity, 
the magnetization, as compared to the technically more 
difficult evaluation of the Kubo formulai^ or the mea- 
surement of time-dependent currents.— Moreover, we can 
control the deviation of the initial state from the ground 
state, thus scanning the regime of systems substantially 
driven out of equilibrium. 

While we show that our results for the spin-1/2 chain 
with an exchange anisotropy in the regime of small per- 
turbations over the ground-state are consistent with the 
picture established by analyzing Drude weights at zero 
temperature, namely ballistic transport in the massless 
and diffusive transport in the massive regime,— we, in 
particular, argue that this also applies to systems far from 
equilibrium. The dynamics is further sensitive to the 
overall filling, or average magnetization, as expected from 
linear-response theory results for the high-temperature 
hmit^iffi^^iS 

Beyond the anisotropic spin-1/2 chain with nearest- 
neighbor interactions only, we further consider noninte- 
grable systems such as the frustrated chain and the two- 
leg spin ladder. As a result, we find that in massive 
phases, the dynamics is typically diffusive, while in the 
massless one of the frustrated chain, the zero tempera- 
ture dynamics are ballistic. 

Transport in the XXZ chain has previously been 
studied in Ref. [13 using tDMRG, there following the 
time-evolution from a highly excited initial state of the 
= I I ... 11 ... I) form. The long-time behavior 
of the magnetization was found to be correlated with 
the phase transition from easy-plane to easy-axis sym- 
metry. Further, the expansion of particle density packets 
of nearly Gaussian form has been looked at with tDMRG 
in the context of ultra-cold atomic gases )^ modeled with 
the ID Bose-Hubbard, as well as for short pieces of in- 
teracting spinless fermionsi^ 

The plan of the paper is the following. In Sec. [Til we 
define the spin models studied and we describe our nu- 
merical method, the tDMRG. We further motivate our 
definition of diffusive transport by discussing the solution 
of the diffusion equation in Sec. IIIII Section IIVI details 
the preparation of initial states. In Sec.|Vl we study the 
magnetization dynamics in the XXZ chain, with ballistic 
transport in the massless regime, and diffusive transport 
in the massive regime. Section IVll summarizes our results 
for two nonintegrable ID systems, the frustrated chain 
and the two-leg ladder. We conclude with a discussion in 
SecEHl 
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II. MODEL AND METHOD 
A. ID Spin-1/2 systems 

Here we will first concern ourselves with the integrable 
XXZ chain: 

L-1 

H = Y.^i^^^tS-+, + H.c.) + ^StSf^^ , (1) 

i=l 

where and// = x,y,z are the components of a spin- 
1/2 operator acting on site i and 5*^^ = S'f ± iSf are 
the lowering/raising operators, respectively. We de- 
note the number of sites by L and we introduce an ex- 
change anisotropy A. Equation ([1]) can be re-expressed 

in terms of spinlcss fermions cf^ through the Jordan- 
Wigner transformation;^ 

L-1 

H = ^[2(4q+i+i/.c.)+AK-l/2)K+i-l/2)] , (2) 

1=1 

with rii = clci. Setting A = results in a nonintcract- 
ing system. If not mentioned otherwise, we impose open 
boundary conditions. We denote the filling factor with 
n. The local magnetization is given by Mi = , and the 
total magnetization is S'' = J2i = ^C*^ ^ 1/2)- 

The ground-state phase diagram of the XXZ chain 
(see, e.g., Ref. [6^ and references therein) exhibits quan- 
tum critical points at A = ±1. A critical phase covers 
the |A| < 1 region, while the ground state for A > 1 
exhibits antiferromagnetic order. The region A < — 1 
has a ferromagnetic ground-state, yet we will restrict the 
discussion to A > 0. The model is integrable through 
the Bcthe-ansatzj^ 

In the second part of this work, we will focus on 
two nonintegrable models with isotropic interactions (i.e., 
A = 1), the frustrated chain and the two-leg ladder (for 
a review on these models, see, e.g., Refs. I66ll68h . Both 
models can be understood as limiting cases of a single 
Hamiltonian that incorporates a dimerization 5 and a 
frustration a: 

H=j [(i+i-imsrsr+i+a^f^r+^i- (3) 



The frustrated chain corresponds to (5 = and a > 0, 
while the two- leg ladder is the 5 = 1 limit. In the latter 
case, we identify the coupling along legs as J\\ = aJ and 
the coupling along rungs as = J . 

The frustrated chain features a quantum phase transi- 
tion at a R:! 0.241, separating a gapless phase from a mas- 
sive onei^ The spectrum of the two-leg ladder is gapped 
for any J^/J\\ > O^i^ 



B. Methods 

For the time-evolution of the noninteracting case (A = 
in Eq. ([2|)), we use exact diagonalization which allows 
us to treat large systems. For all interacting cases, we em- 
ploy the adaptive time-dependent DMRG metho d^^'^^i^° 
with a Krylov-space-based time-evolution scheme22i2i in 
the space of matrix product states.— The control param- 
eters arc the discarded weight 6p and the time-step St. 
Our simulations are canonical ones as we work in sub- 
spaces with a fixed total or particle number, respec- 
tively. We have performed an extended error analysis 
by (i) comparing ED and tDMRG results in the nonin- 
teracting case and (ii) by performing several runs with 
different time-steps and discarded weights at representa- 
tive parameters in the interacting case. Specifically, in 
case (i) , we have analyzed the relative errors in the local 
magnetization 

i=l 

and its variance (see below). It turns out that typically, 
a time step of St = 0.125/J and a discarded weight of 
Sp = 10~^ keeps the relative error SM in M below 10~* 
for a chain of L = 200 sites at half filling and for times 
t < 100/J. 

III. BALLISTIC VS DIFFUSIVE TRANSPORT 

Within linear response theory, one often separates the 
dynamical conductivity cr(w) into a delta-function S{uj) 
at zero frequency and a regular part at frequencies uj > 

Q 59.65 

a{uj) = 2TrDS{uj) + Crrcgular(w) ■ (4) 

These quantities derive from the Kubo formula that 
is based on evaluating current-current correlation 
functions 1^ We repeat that a finite Drude weight in a 
clean, one dimensional system at zero temperature de- 
fines an ideal conductor and thus ballistic behavior, while 
if D = 0, one has an insulator .i^^ The dependence of the 
Drude weight on the exchange anisotropy A in the case 
of the integrable XXZ chain and at zero temperature is 
well-known)^ 



where, in this equation, the anisotropy is parameterized 
through A = cos(i^). We thus have D > for |A| < 1, 
featuring a discontinuous drop to zero at A = 1. Due 
to the excitation gap in the massive regime A > 1, we 
have a true insulator with (Jdc = that can only trans- 
port magnetization once the gap has been exceeded by a 
sufficiently large external perturbation. 
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FIG. 1; (Color online) Density profiles {ni(t — 0)) = Mi{t) + 1/2 in the initial state (time t — 0), induced by (a) a Gaussian 
magnetic field Eq. (|lip {Bo = J, as = 5); (b) a box-shaped magnetic field, Eq. (|12p {Bo = J, sb = 5); and (c) application of 
'^L/2+i' ^^1^ filling n = 0.5. Thin solid lines in (a)-(c): A = 0; thick solid lines in (a): A — 0.5, thin and thick dashed 

lines in (a): averaged density {ni{t = 0)) for A = and 0.5, respectively (see the text in Sec. IIVI for details on the averaging 
and Eq. Results at A = are obtained with ED, all others with DMRG. 



Similarly, the Drude weight vanishes in the massive 
phases of both the spin ladder and the frustrated chain, 
while in the massless regime of the latter model D is finite 
at T = 0^ 

Strictly speaking, on all systems with open boundary 
conditions, the Drude weight vanishes identically. Yet, 
it turns out that the corresponding weight is just shifted 
to small but finite frequencies and thus the system 
is expected to still exhibit ballistic and anomalous trans- 
port properties. 

To justify and motivate our way of analyzing ballistic 
and diffusive transport, let us, for pedagogical reasons, 
consider a ID system that obeys the diffusion equation: 



(6) 



Here, p{x,t) denotes, e.g., a particle density, and V the 
diffusion constant. The Green's function associated with 
this equation is, in a d-dimensional setup, given by 



G{x, X, t) 



1 



{AirVtY/^ 



Therefore, we can calculate the expectation values 



{x){t) = X and (x^) = |fp + 2dT>t 



(7) 



(8) 



and see that the variance cr^ = 
t for normal diffusive transport. 



(x^) — x^ is linear in 
On the contrary, for 
ballistic dynamics, one expects the variance cr^ to grow 
quadratically in time, as is well-known from elementary 
quantum mechanics for free particles. 

Given a distribution of Mi{t) = {jii{t)) — 1/2 at a time 
t, we find it most straightforward to compute the vari- 
ance from the corresponding particle density distribution 
(ni(t)) = Mi{t) + 1/2 as this is a positive quantity. We 
then compute the variance from 



1 



(L/2) 



(9) 



where /i„ is the first moment of the normalized distri- 
bution {rii). Note that we normalize (n^) on the actual 
number of fermions rather than the system size. 



IV. PREPARATION OF INITIAL STATES 

We consider three different ways of preparing the ini- 
tial state, which we now illustrate for the case of the 
XXZ chain, Eq. (IT|), and at half filling n = 0.5, i.e., at 
5^ ~ 0. With the exception of Sec. lVB 2l our simulations 
are always performed in subspaces with these quantum 
numbers. 

Technically, we add a term 



iJs - - ^ BiSi 



(10) 



to the Hamiltonian. By choosing the _B,; appropriately, 
this realizes (i) a Gaussian magnetic field, (ii) a box- 
shape magnetic field, both applied during a ground-state 
DMRG run, and turned off at time i = 0: 



B,{t = Q) 
B,{t = 0) 



Ba exp(-(i-zo)V2<T|) 
L- SB 



e h + 



L + si 



(11) 
,(12) 




FIG. 2: (Color online) Average deviation 5h from half filling 
[see Eq. (|14|l ] in the initial state prepared by applying a Gaus- 
sian magnetic field [Eq. (|11|| ] with crs = 5 as a function of Bo 
for A = 0,0.5, 1 and 1.6. 
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FIG. 3: (a), (b) Time-evolution of the density (magnetization) 
profile (riiit)) — Mi{t) + 1/2 for a Gaussian initial state with 
Bo = J, as = 5 for (a): A = (ED data) and (b): A = 0.5 
(DMRG data, L = 200 sites). 



is the Heaviside function. In Eq. pTjl . ctb is the 
variance of the external Gaussian field, Bq its amplitude, 
and we set = L/2 + 0.5, while in Eq. sb denotes 
the width of a box in which a constant field Bq is applied. 

The third initial state (iii) is realized by finding the 
ground state of a system at filling n = (L/2 — 1)/L and 
then applying a single spin fiip Sf on a site i. The time- 
evolution is then performed at half filling, as in cases (i) 
and (ii). 

The typical shape of the induced density (?ii(t = 0)), 
or equivalently, the magnetization profile is illustrated in 
Figs. [Ua)"(c) at A = (all panels) and A = 0.5 (panel 
(a) only], for the three initial states (i)-(iii), respectively. 
Inherent to the fermionic nature of the model Eq. ([2]), we 
first observe Friedel oscillations with a 2fci? period, where 
kp = 7r/2 is the Fermi momentum at half filling. Sec- 
ond, there are slower spatial oscillations, that are more 
evident in the case of A = 0.5 [Fig. [Ub)]. These oscil- 
lations' characteristic wave-length depends, as we have 
checked, on ub as well as the system size. As we shall 



see later, for the purpose of qualitatively analyzing the 
time-dependence of the variance, the presence of the long- 
ranged oscillations is irrelevant, as the dynamics stems 
from the cental peak dispersing, while away from the 
center, the oscillations contribute subdominantly to the 
time-dependence after turning off Hb ■ Wc will neverthe- 
less sometimes find it illustrating and useful to work with 
a density 



{n,{t)) = [(n2,-i(t)> + ("2.W)]/2 



(13) 



averaged over adjacent sites with i = l,...,L/2 and 
X = 2i — 1/2. To recover the variance of the non-averaged 
density, we multiply by a factor of 2. The aver- 
aged density {hi{t = 0)) is plotted with dashed lines in 
Fig. [IJa) for A = and 0.5, and we see that this averag- 
ing results in quite smooth curves. 

For {hi{t = 0)) and a Gaussian external Bi, we can 
now further address the question whether the induced 
density profile follows a Gaussian as well. We find that 
this is the case at small A, in good approximation. In 
the massive phase A > 1, deviations of {ni{t = 0)) from 
a Gaussian profile are substantial. We will nevertheless 
refer to initial states prepared with Eq. (fTT|l as Gaussian 
initial states throughout. 

As we go from A = into the massive regime A > 1, 
a gap opens, and qualitatively, despite the field Bi being 
inhomogeneous, we expect the existence of the gap to 
affect the average deviation from half filling at a given 
set of {Bq^ob)- To illustrate this point, we display this 
average deviation 5n, defined as 



. Y.^{n.)-l/2Y/N., 



(14) 



in Fig. © (A = 0,0.5,1,1.6). For A = 1.6, we observe 
a steep increase in 5h at Bq/J ~ 0.4, and we will use 
Bq/J> 0.45 for the time-evolution, as below this value, 
little dynamics in the time-evolution is seen. Note that 
Fig. [5] has a logarithmic scale on the dn axis, and below 
Bq/J « 0.4, Sh oc exp(c/i3o) in the case of A = 1.6. We 
find c = 0.4J and this roughly coincides with twice the 
spin gap for this value of 



V. THE XXZ CHAIN 

After detailing the way of preparing initial states, we 
now come to the analysis of Mi(t) = {ni{t)) — 1/2, fo- 
cusing on the variance a\j{t). In this section, we will 
first discuss cr\[{t) in the massless regime in Sec. IV Al 
and show that we clearly see ballistic behavior with 
c\[{t) = const -f D2t^ ^ independently of the initial state. 
We will then analyze the dependence of the coefficient 
0-2 on Bf), scanning the full range of perturbations from 
a linear one with D2 oc B^ to the largest perturbations 
possible. These regimes are distinct by a different finite- 
size scaling behavior, to be discussed below. As we illus- 
trate in the case of a Gaussian magnetic field Eq. pip . 
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FIG. 4: (Color online) Time-evolution of the density {ni{t)} 
(or magnetization profile Mi(t) = {ni(t)) — 1/2) for a Gaussian 
initial state with _Bo = J, era = 5 at A = 0.5 (DMRG data, 
L = 200 sites): Snap-shots from Fig.[3]at times t J = 0, 15, 25. 

ballistic dynamics is found for A < 1. This holds inde- 
pendently of the actual choice of (Bo,aB) and thus also 
far from equilibrium, except for the most extreme initial 
states considered here (see the discussion in Sec. IV A 21 
below) . 

We then, in Sec. IV Bl discuss the transition from bal- 
listic to diffusive behavior which is expected to occur at 
A = 1. Our data are consistent with this picture, as we 
find strong evidence for diffusive transport for A > 1.5. 
We will present results for several So at A = 1.5 to sub- 
stantiate that the observation of diffusive transport is 
independent of the initial state. 

A. XXZ chain: Massless regime 

1. Time- dependence of the variance 

We now turn to the analysis of Mi{t) = {ni{t)) — 1/2 
in the massless regime. Figures El^a) and (b) show 
(ni(t)) = Mi{t) + 1/2 for A = and A ^ 0.5, respec- 
tively. The initial density profile first melts and then, at 
times t J ^ 5, splits into two packets that travel into op- 
posite directions!^ Figure [5] shows snapshots of {ni{t)) 
at times < J = 0, 15, 25 for A = 0.5. It is noteworthy 
that, while substantial oscillations are present far away 
from the central peak, these oscillations are frozen in and 
do not contribute to the increase in a^j since, far away 
from the center of the chain, H ^ H + Hb- 

The variance (j\[{t) — (j\j{t = 0), plotted vs. time, is 
displayed in Fig. EJa), for A = (solid line), A = 0.5 
(dashed line), and A = 1 (dotted Hue) for the evolution 
from an initial state of the type (i), enforced by a Gaus- 
sian magnetic field with Bq = J/2 and gb = 5. The 
circles represent the time-evolution of a\[ (t) — [t = 0) 
computed from the averaged density {n{t)) (sec Sec. IIV[) 



for A = 0, 0.5, 1. For A = and 0.5, the averaged results 
very well coincide with a\,j{t). For the purpose of char- 
acterizing ballistic or diffusive behavior, it therefore does 
not matter whether the pure fermionic density {ni(t)) 
or the averaged quantity {ni(t)) is used. In what fol- 
lows, we will present results extracted from the former, 
unless stated otherwise. We mention, though, that the 
quantitative difference in the variance extracted from the 
averaged as compared to the bare density becomes more 
pronounced the larger A and the smaller Bq is. This be- 
comes evident in the case of A = 1, included in Fig.[5ja). 

The key observation from Figs. EUa) and (b) is the 
quadratic increase of the variance with time observed for 
A = 0,0.5 and 1, which confirms the expected ballistic 
behavior in the critical regime. For the isotropic chain 
(A = 1), we find that the best fit of a power-law to 
<it) 7 <it ~ 0) yields aj.it) - ^^.(t = 0) « 
which is thus slightly below the behavior expected for 
ballistic transport. However, a deviation of just one per- 
cent from the expected exponent of 2 is very much within 
the accuracy of our numerical calculations. 

Moreover, this behavior, as we show in Fig. EJb) for 
the example of A = 0.5, is independent of the shape of 
the initial state: all three types of states studied here - 
Gaussian field, box shape, and application of - result 
in ballistic dynamics at half filling. 



2. From small to large perturbations 

So far, we have worked at a fixed pair of Bq and cfb- 
We will now explore the dependence of the dynamics on 
how far the initial state is perturbed away from the actual 
ground state. This deviation is measured through the av- 
erage deviation Sh from half filling (compare Fig.[21). For 
that purpose, we fix ctb = 5, and perform a set of simula- 
tions at different Bq, at A = 0,0.5 and 1. Independently 
of Bo, we always find a quadratic increase of cr^^(i) with 
time, i.e., cr]^{t) = const -I- D2{Bo)t'^. As a main result, 
we therefore conclude that the qualitative behavior of the 
dynamics is independent of the external perturbation, it 
is always ballistic in the massless regime. 

On a more quantitative level, it is instructive to plot 
D2 vs Bq, shown in Fig. [D^c) for A = and 0.5, as 
it allows us to distinguish two regimes: first, the linear 
regime, in which D2 is linear in Bq with D2 L/2 ^ f{L) 
at a fixed Bq. For A = and 0.5, the linear regime 
extends up to i?o/J ~ 0.5. Second, at larger Bq effects 
of both the band curvature and the finite band-width 
start to play a role, with significant finite-size effects, as 
illustrated in Fig. [6] for Bq — 2 J. In the case of A ~ 
and Bq = 2 J, we are able to access system sizes of i « 
10^, and at large L, the scaling is of the form D2 L/2 oc 
1/L, allowing for an extrapolation to i — > 00. In the 
interacting case, the accessible system sizes are too small 
to establish such scaling and we thus have not attempted 
any extrapolation in the case of A = 0.5. 

Let us next discuss the limiting cases of first, the linear 
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FIG. 5: (Color online) (a) Time-dependence of the variance aM{t) — (j\i{t = 0) for A = (solid lines) and A = 0.5 (dashed 
lines), both at Bo = J/2 and gb ~ 5. Circles denote the variance o"m {t)-alj(t = 0) of the averaged density {hi{t)), Eq. ([13}. 
(b) Time-dependence of the variance crf/(t) — (J%j(t — 0) for the evolution from a box-like initial state (solid line), after the 
application of S^^^j (dot-dashed line), and from a Gaussian-state (dashed line) (all at A = 0.5). The squares denote the perfect 
fit of all sets to <jfj{t) — cr\j{t = 0) = D-^t^ . (c) _Bo-dependence of the coefficient D2 of the time-dependent variance for 
A = (solid lines), A = 0.5 (dashed lines), and A = 1 (circles) at a fixed as = 5 and L = 200. The dotted lines display 
D2 L/2 — 4:BoJ. Inset: D2{Bo) in the limit of large Bo = 10' J for A = 0. The thin dashed line is the result for the expansion 
from a Fock state of width L/2 with (rii) = 1 at the center of the chain. 



regime, i.e., D2 oc Bq, and second, the limit of Bq — > 00. 
Starting with the former, the linear regime, D2L/2 = 
•yBQ, we find that the prefactor is 7 « 4J for values of 
A < 0.5. At large A, this reduces to 7 « 3.4J as the 
results for A = 1 displayed in Fig. [5Uc) show (circles). 
The interpretation of D2 being linear in Bq for A < 1 is 
based on the observation that the area ^pcak(-Bo) under 
the initial Gaussian-like magnetization profile increases 
linearly in Bq , which wc can strictly confirm in the A = 
case and chains with up to L = 1000 sites. This analysis 
requires an estimate of the background density, which 
can the best be done in the A = case, but suffers 
from finite-size effects at a nonzero A. In the A = 
case, wc find that 7-Bo/Apeak(^o) = 4J^ and for A = 
0.5, 7i?oMpcak(-So) ~ 1.7J^. We may therefore conclude 
that in the linear regime, j Bq = Apesk{Bo)v^{A) , where 
Vg{A) is the group velocity (see, e.g., Ref . 1781) : 



TT sin(zy) 



(15) 



with A = cos(:^). From our numerical data for L = 200 
sites and at A > 0.5, we obtain Z?2 ^/(SApoak) > 'Vg{A) 
due to finite-size effects, but qualitatively, D2 L/ {2 Apeak) 
increases with A at a fixed Bq in the linear regime as 




FIG. 6: (Color online) Finite-size scaling of D2 for A = and 
0.5 in the large Bo regime (Bo = 2J, as = 5). 



expected from Eq (15) for the behavior of the velocity Vg. 
At sufficiently large Bq , we expect the band-curvature to 
play a role as well which we are able to verify in the case 
of A = 0. D2 L/(2^poak) then decreases below Vg{A = 0) 
as Bq increases. 

As for the limit of large Bq ^ J, note that since Mi < 
1/2 and since we work at fixed half filling of the full 
system, the most extreme initial state that, on a fixed 
system size L, Bq ^ 00 drives the system into is a Fock 
state I/) with {ni{t = 0)) = 1 for L/4 < i < 3L/4 and 
zero otherwise. Wc thus follow the evolution from such 
a state as well as the evolution from states with very 
large Bq ~ 10' J (j = 1, . . . , 30). The inset in Fig. E^c) 
shows that as Bq increases, D2{Bq) indeed approaches 
the value found for the limiting state |/) for Bq > 10^^ J 
in the example of A = 0. In that limit, D2 = 1/2 J^. 

While for the parameters of the main panel of Fig.OJc), 
i.e., A ~ 0,0.5,1 and Bq/J < 2.4, the variance always 
follows a power law with exponent two, curiously, this is 
not the case for the aforementioned Fock states encoun- 
tered in the i?o — *■ 00 limit and A = 0.5,1. There, we 
find exponents that are consistently below two for A > 
and L = 40, 80, 200. This behavior was also observed in 
Ref. [53 for the evolution from similar Fock states. A full 
analysis of the evolution from Fock states will be pre- 
sented elsewhere. 

We mention that the time-dependent evolution from 
Fock states or, more generally, the evolution of particles 
originally trapped in a confined region into an empty lat- 
tice, has been intensely studied with numerical methods 
in Refs.lliliaEllIll, for the cases of hardcore boson s ' , 
soft-core bosonS)^ and the interacting two-component 
Fermi gasi^ These studies have a fundamental inter- 
est in out-of-equilibrium phenomena, with a perspective 
onto experiments with ultra-cold atoms in optical lat- 
tices. Among these we mention the experimental efforts 
directed at detecting Anderson localization in cold atom 
gases, precisely by utilizing such expansion sct-ups42i'^ 
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FIG. 7: (Color online) (a) Transition from ballistic to diffusive transport: — cr^/(f — 0) for A = 1 (thick solid line) 

A = 1.2, 1.3, 1.4, 1.5 (thin solid lines). The parameters are Bo = 2J, as = 5, L — 200. The thick dashed line shows the 
asymptotic behavior for A = 1.5, i.e., aM{t) = const + Dit. (b) Variance for A = 1.5 and several Bo/ J — 1, 1.5, 2 (L — 200, 
ub ~ 5). (c) Variance for A = 1.5 at = 0, 5, 10, 15, 20 (see the legend) and -Bo = 2J, as = 5: ballistic transport is restored 
away from half-filling. 



Finally, while all results discussed here were obtained 
from chains with L = 200, we stress that we have care- 
fully studied the finite-size scaling of Z?2(see Figs. [5jc) 
ancini). By plotting D2 L/2 in Fig.EJc), wc account for a 
trivial sizc-dcpcndcnce. and curves obtained from differ- 
ent L but the same A indeed collapse onto a single one 
in the linear regime. At larger Bq, D2 tends to decrease 
with system size L as less particles per Bq can be accu- 
mulated in the initial Gaussian peak, due to Sf < 1/2 
(see Fig. |6]). 

B. XXZ chain: Massive phase 

1. Transition from ballistic to diffusive behavior 

Linear-response theory predicts a sharp transition from 
ballistic spin transport, characterized by a finite Drude 
weight D, to diffusive spin transport at A = 1.— 
We have carried out several simulations with different 
A = 1, 1.2, 1.3, 1.4, 1.5 at fixed parameters Bq = 2J and 
ub = 5 to check whether an analysis of the variance cap- 
tures this transition. As the data for <j\i{t) displayed in 
Fig. Elja) show, we clearly find a linear increase of (j\i{t) 
at large times in the case of A = 1.5, which, in the sense 
of Sec. mil we interpret as evidence for diffusive trans- 
port. At smaller 1 < A < 1.4, our data do not allow for 
this conclusion, yet at least we can state that the data 
for A = 1.3 and 1.4 do not follow a power-law, indica- 
tive of non-ballistic transport. Also note that there is 
no contradiction with Eq. ([5]) as it is is well-known that 
finite-size effects in the vicinity of A = 1 are severe and 
come along with a logarithmically slow convergence with 
system size of quantities such as the Drude weight or the 
spin stiffness jiSii^ We suspect that larger system sizes, 
hence access to longer simulation times, are necessary to 
fully capture the sharp transition from ballistic to diffu- 
sive behavior at A = 1. We stress that with iSg = 2 J, 
we work with highly perturbed initial states, and thus 
the observation of diffusive transport for A > 1.5 is a 
nontrivial one, going beyond the case studied in Ref. [sol . 



Our results for A > 1.5 thus establish an example of dif- 
fusive dynamics with a\j {i) = const -I- Di t in this model 
for the out-of-equilibrium situation. Note that a recent 
tDMRG work on transport in spin chains incorporating 
baths has reported similar results, derived from current 
and spin profiles in the steady statei^ 

It is noteworthy that at short times, obviously, the 
dynamics is always ballistic, independently of A, as can 
be seen in Fig.[71[a). Even the value of (j\[{t) is roughly 
the same for all A at short times. In the long time limit, 
which is the relevant one to characterize the system as 
diffusive or ballistic, we find that cr\,j{t) systematically 
decreases with increasing A. The reason is that, in the 
A/J cxD limit, no dynamics is possible at all. 

We have further studied the dependence on Bq at IS. ~ 
1.5. In Sec. IIVI we hinted at the fact that at A > 1, a 
reasonably large Bq is necessary to observe a significant 
change in the magnetization profile over the time scales 
simulated. We attribute this to the existence of a spin 
gap in the antiferromagnetic phase A > 1 and focus on 
-Bo ^ J ■ Our results of several runs at A = 1.5, scanning 
the Bf) dependence, are displayed in Fig. [7{b). We find 
that the variance increases with Bq, which we mainly 
attribute to more particles accumulated in the central 
peak. Furthermore, the figure suggests that the time- 
scale at which diffusive behavior sets in depends on Bq 
such that for larger Bq, we are able to observe a linear 
increase of cf\i{t) earlier in time. 



2. Restoring ballistic transport 

While so far we have restricted ourselves to the case 
of half filling, we now address the magnetization dynam- 
ics at incommensurate filling. The initial states are now 
created by applying the external fields in subspaces that 
already have a finite magnetization. , i.e. 7^ 0. Results 
for the variance at A = 1.5 are displayed in Fig.[7Jc). The 
z-component of the total spin is 5*^ = 0, 5, 10, 15, 20. We 
find that any nonzero is sufficient to render the dy- 
namics ballistic again, on the time scales accessible to our 
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FIG. 8: Two-leg spin ladder with = J± [a = 1, S — 1 in 
Eq. dS])]: Contour plot of {ni{t)) = Mi{t) + 1/2. In the initial 
state, a Gaussian external field Bi with Bo = J and as = 5 
is applied [see Eq. ([11])]. DMRG data, L = 200 sites. 

simulations. According to our data, the variance follows 
alj{t)-alj{t = 0) cx with a « 2.05±0.02. This obser- 
vation is in agreement with the infinite-temperature be- 
havior of the Drude weight )^^i^° which in the XX Z chain 
is finite at any A away from half filling. Note that A > 1 
and 5^ > but below saturation is in the easy-plane 
phase of the XX Z chain, with gapless excitationsi^ Note 
that the prefactor in cr|f (t) - (T^(t = 0) = D2(5"^)r de- 
pends on 5*^ in a non-monotonous way: it is the largest 
at around S*^ = 10 and then decreases as saturation is 
reached. 



VI. NONINTEGRABLE MODELS 

We finally move on to the discussion of the dynam- 
ics in two nonintegrable models, the two- leg ladder and 
the frustrated chain, both limiting cases of Eq. ([3]). Nu- 
merical studies of the /ii(7/i-tempcrature limit, based on 
the Kubo formula, conclude that spin and thermal trans- 
port in the massive phases of these models (see Sec. HI]) 
are normal, with a vanishing Drude weight i^ii^i^ii^ The 
conclusions on the massless phase of the frustrated chain 
are not unambiguous r^"^*^ and it has been pointed 
out that the energy-current operator, to first order in the 
next-nearest-neighbor interaction aJ, is conservedi^ 

One scenario is that the high-temperature Drude 
weight vanisheS )^°i^^ while it is still possible to find 
anomalous transport properties in the low-temperature 
regime, e.g., in the form of a peculiar low-frequency be- 
havior of crroguiar(^)- Exact diagonalization results show 
that the Drude weight is finite at zero temperature in the 
massless phase of the frustrated chain ^iiS 

We here use the approach outlined in the previous sec- 
tions to show that the zero-temperature dynamics of two- 
leg ladders and frustrated chains with a spin gap is of 



diffusive nature. To this end, we prepare initial states 
with Gaussian magnetic fields Eq. pT]) . We emphasize 
that, in the case of the spin ladder, both sites on a rung 
experience the same field. In these two cases, and similar 
to the discussion of initial states for A > 1 in the XX Z 
chain (compare Fig. [2]m Sec. lIVp . the amplitude of the 
Gaussian field, Bq, needs to be large enough to induce 
a substantial perturbation in the magnetization Ali that 
will actually propagate through the system. We thus here 
probe the magnetization dynamics and transport at large 
external perturbations. 

Starting with the example of a spin ladder with Jx = 
J|| , we display the magnetization profile Mi{t) = {ni{t)) — 
1/2 in Fig. [8] as a contour plot. The time-dependence of 
the corresponding variance is shown in Fig. OJa) [solid 
line, squares], and we find a linear increase in cr^/(i) for 
times t >17/ J, clearly establishing the notion of diffusive 
dynamics in the ladder system. 

A more involved picture emerges in the case of the 
frustrated chain. For this model we present results for 
a = 0.2 (circles) and a = 0.4 (stars) in Fig.lHKb). While 
on the time scales simulated, the variance for the a = 
0.2 curve perfectly follows the form (T^/(t) = const + 
D2t^ (a least-square fit to this function is displayed by a 
thin solid line), in the case of a = 0.4, we observe that 
the data do not follow a power law cr^/(i) = const + 
Dajf', which supports the notion of a time-dependent 
crossover from ballistic to diffusive dynamics. In fact, the 
numerical results yield a variance that clearly increases 
linearly in time for t > 30/ J. Note that the transition 
from the massless to the massive phase in this model is of 
the Beresinski-Kosterlitz-Thouless type^^i^i^^ with an 
exponentially growing correlations lengths as the critical 
point acrit is approached from a > ctcrit- This renders 
it very difficult to see a sharp transition in the transport 
behavior using exact diagonalization or DMRG as acrit 
is crossed. 

We mention, though, that for other values of a > 
acnt ~ 0.241 (results not shown here), on similar time- 
scales, no diffusive behavior is seen. Moreover, the 
time-scale at which diffusive dynamics emerges seems to 
strongly depend on Bq, i.e., on how far the initial state is 
perturbed over the actual ground state. The qualitative 
trend is that the larger Bq, the faster diffusive transport 
is established. Unfortunately, the larger Bq, the worse is 
the entanglement growth which renders tDMRG simula- 
tions more difficult (sec, e.g., Ref. i87i). 

Keeping in mind these remarks, we are in a posi- 
tion to conjecture that in general, in massive phases 
of one-dimensional spin models, spin transport is diffu- 
sive. Combined with the existing results for the high- 
temperature limit our work suggests that this ob- 
servation applies independently of temperature. 

Note that a recent exact diagonalization study^ has 
promoted a different behavior, namely evidence for bal- 
listic spin transport at zero temperature in the frus- 
trated spin chain at a = 1. This conclusion is based 
on the presence of certain oscillations (dubbed a bounc- 
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FIG. 9: (Color online) Time dependence of the variance 
o"m(*) — o"Af = 0). (a) Spin ladder, solid line with squares, 
compiled from the data of Fig. |S] Bo = 1.7J. (b) Frustrated 
chain with a = 0.2 (circles, Bo = J) and a frustrated chain 
with Q = 0.4 (stars, Bq = 1.5J). In all cases, ob = 5. The 
dashed lines in (a) and (b) are a fit of const+Di t to the vari- 
ance of the spin-ladder and the frustrated chain with a = 0.4 
at long times, while the thin, solid line in (b) is a least-square 
fit of D2 to the result for the frustrated chain with a = 0.2. 
DMRG data, L = 200 sites. 



ing hehavior^^ in Miit) in the evolution from an initial 
state with all spins pointing up (down) in the left(right) 
part of an open system (compare Refs. [27ll57h . We beheve 
that our analysis of the variance is more quantitative, 
which may imply that the oscillations in Miit) reported 
on in Rcf. [s^ are possibly of a different origin. A recent 
study of quantum quenches in the XX Z chain proposes 
that oscillations seen in the order parameter are related 
to the quantum phase transition at A = 1 (Ref. [s^, see 
also Ref. EF 



VII. SUMMARY AND DISCUSSION 

In this work we studied the nonequilibrium magneti- 
zation dynamics in one-dimensional spin models at zero 
temperature using the adaptive time-dependent DMRG 
method on system sizes as large as L = 200 sites. We 



considered several models: the integrable spin-1/2 XX Z 
chain, the frustrated chain, and the two-leg spin ladder. 
Based on the analysis of the time-dependence of the spa- 
tial variance of the magnetization during the time evolu- 
tion starting from initial states with an inhomogcneous 
magnetization profile, we conclude that in the critical 
regime of the XX Z chain, the magnetization dynamics is 
ballistic. In contrast to that, in the massive regime, our 
results indicate diffusive transport at half filling, while 
ballistic transport is restored away from half filling. A 
major aspect of our work is that we scanned the entire 
regime going from small to very strong perturbations over 
the ground state. This substantially extends previous 
studies of linear-response functions as we clearly enter 
into a regime with the system driven out of equilibrium. 
In the case of the masslcss regime of the XXZ chain, 
ballistic transport is seen for substantially perturbed ini- 
tial state, while for the most extreme initial states, i.e., 
pure Fock states, we still find a power-law for the time- 
dependence of the variance, but, on the times scales sim- 
ulated, with an exponent below twoi^ 

As for the nonintcgrable models, the frustrated chain 
and the ladder, the numerical data clearly support the 
notion of diffusive dynamics in the ladder system. In 
the case of the frustrated chain, our data are consistent 
with a transition from ballistic to diffusive behavior as 
the quantum critical point acHt ~ 0.241 is crossed. In 
the limit of small perturbations, this result confirms with 
the general picture that masslcss phases support ballis- 
tic, and massive ones diffusive dynamics, at zero temper- 
ature and irrespective of integrability4 Overall, a differ- 
ence between the low- and high-temperature behavior is 
then evident: numerical results for the high-temperature 
limit^ consistently support the notion of vanishing Drude 
weights in nonintcgrable models, and thus normal trans- 
port behavior, irrespective of what the ground state 
phases are. Conversely, exact diagonalization studies find 
a finite Drude weight in the gapless phase of the frus- 
trated spin-1/2 chain at zero temperaturei2iiSi2£ In the 
low, or more extremely, zero temperature case, effective 
low-energy theories are expected to give a valid descrip- 
tion, which typically predict diverging transport coeffi- 
cients of clean spin systems (see, e.g., Ref. [911) . As for 
the heat transport measurements on spin chain and lad- 
ders experiments (see Refs. [siUb^ for a survey), a domi- 
nant magnetic contribution is usually evident in the high- 
temperature regime, where the validity of effective low- 
energy theories for the description of transport is not 
obvious. 

Finally, the approach of distinguishing ballistic from 
diffusive transport by analyzing the spatial variance of 
a density like quantity could be instrumental in charac- 
terizing ultracold atomic gases in optical lattices as well. 
There, one typically realizes the expansion of particles 
into an empty lattice, and experimentally it is possible 
to measure the expanding cloud's radius. It would thus 
be very interesting to identify conditions for ballistic as 
compared to diffusive dynamics for model systems typi- 
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cally encountered in ultracold atomic gases such as the acknowledge support from the Deutsche Forschungsge- 
Bose-Hubbard model or a two-component Fermi gas in meinschaft through FOR 912. 
an optical lattice. 
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